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Abstract. Effect of stellar electromagnetic radiation on motion of spherical dust 
^ ' particle in mean-motion orbital resonances with a planet is investigated. Planar 

-s_i ' circular restricted three-body problem with the Poynting-Robertson (P-R) effect 

\l ' yields monotonous secular evolution of eccentricity when the particle is trapped 

^^ in the resonance. EUiptically restricted three-body problem with the P-R effect 

("^ ' enables nonmonotonous secular evolution of eccentricity and the evolution of ec- 

^\ ^ centricity is qualitatively consistent with the published results for the complicated 

case of interaction of electromagnetic radiation with nonspherical dust grain. Thus, 



it is sufficient to allow either nonzero eccentricity of the planet or nonsphericity 
of the grain and the orbital evolutions in the resonances are qualitatively equal 
for the two cases. This holds both for exterior and interior mean-motion orbital 
resonances. Evolutions of longitude of pericenter in the planar circular and ellipti- 
cal restricted three-body problems are shown. Our numerical integrations suggest 
that any analytic expression for secular time derivative of the particle's longitude 
of pericenter does not exist, if a dependence on semi-major axis, eccentricity and 
longitude of pericenter is considered (the P-R effect and mean-motion resonance 
with the planet in circular orbit is taken into account). 

Change of optical properties of the spherical grain with the heliocentric dis- 
tance is also considered. The change of the optical properties: i) does not have any 
significant infiuence on secular evolution of eccentricity, ii) causes that the shift 
of pericenter is mainly in the same direction/orientation as the particle motion 
around the Sun. The statements hold both for circular and noncircular planetary 
orbits. 
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1. Introduction 

Orbital motion of dust particles in the zones of mean-motion orbital resonances with a 
planet is intensively discussed since the paper by Jackson and Zook (1989) was published. 
A body is in resonance with a planet when the ratio of their mean motions (mean motion 
n = Stt/T, where T is orbital period) is the ratio of two small integers. Besides important 
gravitational attraction of the Sun and the planet, moving usually in circular orbit, also 
the effect of solar electromagnetic radiation on the particle is considered. Standardly, 
dust particle is considered to be spherically symmetric and, correspondingly, the effect 
of solar electromagnetic radiation is considered in the form of the Poynting-Robertson 
(P-R) effect: e. g., Jackson and Zook (1989), Sidlichovsky and Nesvorny (1994), Beauge 
and Ferraz-Mello (1994), Marzari and Vanzani (1994), Liou and Zook (1995), Liou et al. 
(1995), Liou and Zook (1997). Observations confirming the existence of dust ring around 
the Sun in resonant lock with the Earth are discussed in Brownlee (1994), Dermott et 
al. (1994), Reach et al. (1995). The paper by Liou and Zook (1997) presents secular 
orbital evolution of eccentricity when the particle is trapped in mean-motion orbital 
resonances (Liou and Zook 1997 - Eq. 26; see also Klacka and Kocifaj 2006a - Sec. 9, 
mainly Eq. 83). As a consequence, P-R effect and circular restricted three-body problem 
yield monotonous secular evolution of eccentricity of the particle trapped in a mean- 
motion orbital resonance with a planet. However, qualitatively new result was obtained 
by Klacka et al. (2005), when a general interaction of electromagnetic radiation with dust 
grain was taken into account. Really, arbitrarily shaped dust grain and its interaction with 
solar electromagnetic radiation does not yield, in general, monotonous secular evolution 
of eccentricity of the particle trapped in mean-motion resonances (Klacka et al. 2005a, 
2005b; Klacka and Kocifaj 2006a, 2006b). Thus, an interesting question arises: Does exist 
any other generalization of the standardly used access (P-R effect + planet in circular 
orbit) which can produce results qualitatively consistent with nonspherical grains? 

We concentrate on orbital evolution of spherical dust particles in mean-motion orbital 
resonances, in this paper. We take into account nonzero planetary eccentricity, also. 
Besides constant optical properties of the particle, also dependence of optical properties 
on heliocentric distance is considered, in accordance with Kocifaj et al. (2006). Influence 
of this effect is usually ignored in literature. In this paper we concentrate on orbital 
evolution of dust grains inside resonances: the initial conditions are given in the form that 
the particle starts its motion inside the resonances. Results of our study can be applied 
also for the case of orbital evolution of particles captured in the Earth resonant ring, which 
is well observed (Brownlee 1994, Dermott et al. 1994, Reach et al. 1995). We compare 
our results with the results obtained for nonspherical particles which were published 
elsewhere. It is shown that nonzero eccentricity of the planet can mimic behaviour of 
nonspherical dust grains for circular planetary orbit. 
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Attempts in dealing with restricted elliptic three-body problem can be found also in 
Gonczi et al. (1983). However, the authors ignore (see their Eq. 1) the inertial gravita- 
tional term and the nongravitational term represents only a partial component generated 
by the action of solar electromagnetic radiation (they would not receive, e. g., results on 
orbital evolution presented by Robertson 1937 and Wyatt and Whipple 1950). Thus, 
our equation of motion will be more general, and it will really correspond to physics. 
Moreover, we are motivated by the results obtained during the last ten years, which were 
not known at the time of Gonczi et al. (1983). 

2. Mean-motion resonances with planets 

According to the third Kepler's law we have 

a^ fi^ = G Mq (1 - /3) , 
a%nl = G {Mq + mp) , (1) 

where a, ap are semi-major axes of a particle (characterized by optical parameter /?) and 
a planet with mass mp, n and np are mean motions characterizing revolutions around 
the Sun of the mass Mq. The first part of Eq. (1) uses the fact that central Keplerian 
acceleration is given by the sum of solar gravitational acceleration and radial component 
of radiation pressure acceleration. Eq. (1) yields 

-1/3 



-(!-«'« (^)"0-S) 



ap . (2) 



A particle is in a mean-motion resonance with a planet when the ratio of their mean 
motions is the ratio of two small integers. If the dust particle is in the resonance with the 
planet, we can define the g— th order exterior resonance by the relation np/n = ip + q)/p, 
where p and q are integer numbers. Similarily, the q—th order interior resonance is defined 
by the relation np/n = p/{p + q), where p and q are integer numbers. In terms of orbital 
periods: T/Tp = {p + q)/p for exterior , T/Tp = p/{p + q) for interior resonance. On the 
basis of these definitions and Eq. (2), we can immediately write 

/ I \ 2/3 / \ -1/3 

for the semi-major axis of the dust particle in the q—th order exterior resonance with 
the planet of mass mp. Similar relation can be obtained for the interior resonance. 

3. Equation of motion for spherical dust particle 

Let us consider a spherical dust grain under action of gravitational forces generated by 
the Sun and a planet moving around the Sun. Let the grain is under action of solar 
electromagnetic radiation, too. The considered electromagnetic radiation effect is known 
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Fig. 1. Radiation pressure efRciency factor Q' as a function of heliocentric distance for 
compact carbonaceous dust particles with radii R = 2 fiia and R = 5 /inr. 



as the Poynting- Robertson effect. Equation of motion of the particle is considered in the 
form 



d V GMq (1-/3) ^ GMq 

- — Er ~ P 



V -Er V 

eR + - 

c c 



P 



Gmp 



Lp,TrR^ 



rp 



\r — rp\ 



rp 



i'pr 



7.6 X 10" 



-, 7ri?2[m2] 



(4) 



Att GMqTti c ^" '^P'' mlkg] ' 

where r and rp are position vectors of the particle and the planet (mass mp) with 
respect to the Sun (mass Mq), r ~ \r\, br = r/r, v = dr/dt is velocity of the particle, 
G is the universal gravitational constant, c is speed of light, m is mass of the particle 
[m — AirR^p/i we consider homogenous particles with density p), i? is radius of the 
particle, Q' is the spectrally averaged efficiency factor for radiation pressure and Lq 
is the solar luminosity. As for the Poynting-Robertson effect, we use equation of motion 
derived by Robertson (1937), see also Klacka (1992, 2000, 2004). 

In this paper we consider Eq. (4) as the equation of motion in the elliptically re- 
stricted three-body problem together with the P-R effect. To obtain planetary position 
in elliptical orbit we solve Kepler equation. Initial conditions of the planet and the particle 
always correspond to prograde (counter-clockwise) motion, in our numerical integrations. 
Moreover, in some cases, we admit that optical properties of the particle may change with 
heliocentric distance due to the change of temperature. Change of optical properties is 
described by change of spectrally averaged efficiency factor for radiation pressure Q'. 
We have calculated the value of Q' for particles with radii 2 /im and 5 /im in the range 
of heliocentric distances from to 10 AU. Relation of Q' on heliocentric distance is 
shown in Fig. 1. Method that we used for calculating spectrally averaged efficiency factor 
is taken from Kocifaj et al. (2006) and the results of Klacka et al. (2007). 
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4. Secular evolution of eccentricity in mean-motion resonances with planet 
in circular orbit 

Secular evolution of semi-major axis is characterized by its constant value when the 
particle is in a resonance with a planet. What can be said about secular evolution of the 
eccentricity of the particle in such a situation? In this derivation we will suppose that 
the planet is moving in a circular orbit around the Sun. In this case, we have a special 
gravitational problem of three bodies and small nongravitational forces are also present. 
The gravitational problem is know as the circular restricted three-body problem. At the 
end of the 19-th century F. F. Tisserand found a quantity, which does not change during 
motion of the third body, whose mass is negligible in comparison to the masses of planet 
and the Sun (see, e. g. Brouwer and Clcmcnce 1961). For our case, we can write the 
Tisserand's parameter in the form 



1-/3 n /(l-/3)a(l-e2) 

Ct = + 2 J -i -^-V cos / , (5) 

a \j aj, 

where e is eccentricity of the particle characterized by constant parameter /? and / is the 
inclination of the particle's orbital plane with reference to the plane of the planet's orbit. 
We have to stress that Tisserand's quantity Ct does not change only for the special case 
of the circular restricted problem of three bodies. However, Eq. (5) enables to find secular 
change of the eccentricity of the particle captured in a resonance. For this purpose, we 
will consider / = in Eq. (5). We can write 

dCr _ dCr fda\ dCr fde\ 

for the total time derivative of the Tisserand quantity C't defined by Eq. (5). However, 
according to Eq. (4), time derivatives of semi-major axis and eccentricity of the particle 
are caused by gravitational perturbations of the planet (these terms will be denoted by 
the subscript G) and nongravitational perturbations caused by Poynting- Robertson effect 
(these terms will be denoted by the subscript NG): 
da\ / da \ /da 

^) total " \^)g ^ V^ 

de\ f de\ f de 



"'^ / total \"-^ / G V "'• / NG 

On the basis of Eqs. (6)- (7) we can write 
dCx dCr ( f da\ (da 



dt da \\dt J Q \dt / j^q 

According to Tisserand, gravitational terms alone do not change the value of Ct- 

dC^fda\ aCV fde\ ^ 

da \dtJc ^ de \dtj^ ^' 
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Putting Eq. (9) into Eq. (8): 

d^T^dCr fda\ ^ dCr f de\ 

dt da \dt J jyg, de \dt J ^g. 

If we are interested in secular changes of orbital elements a and e, then the particle's 
stay in the resonance is characterized by the relation (da/dt) = 0: for a function a of the 
property a{T) = a(0) the relation (da/dt) = [l/T) j'^ [da/dt) dt = (a{T) - a{0))/T = 
holds. After averaging over period of the resonant oscillation of semi-major axis Eq. (6) 
reduces to 

d^\^dCT /de\ 
dt / de \dt/ ' ^ ' 

Averaged Eq. (10) and Eq. (11) finally give for the total secular change of the eccentricity 
of the particle 

de\ _ /^\ , dCx/da / da\ 

dt/'\dt/^^ dCr/de \dt/^^- ^ ' 

Calculating partial derivatives of Ct defined by Eq. (5), and using relations for secular 
changes of the semi- major axis and eccentricity for the P-R effect (with assuming constant 
optical properties of particle) (Robertson 1937; Wyatt and Whipple 1950; Klacka 2004) 
/da\ __ GMq 2 + 3e2 
\dt/^a~ ^ c a(l-e2)3/2' ^1^) 

^\ ^_5^Mo t (14) 

di/jvG 2 ^ c a2(l-e2)i/2 ' ^ > 

we get for the total secular change of the eccentricity (under the assumption that Eqs. 
(13)-(14) hold also for the period of resonant oscillation of the semi-major axis) 

dt / c a^ e 

ia/a,f^ (1 - e^f^ ' ^''^ 

see Liou and Zook (1997), Klacka and Kocifaj (2006a, 2006b). 

Eq. (15) determines secular evolution of eccentricity of the spherical particle char- 
acterized by constant values (3 and Q' , if the particle is captured into a mean- motion 
resonance with planet moving in circular orbit. This equation enables to find detail evo- 
lution. If we take some special mean-motion resonance, we already know the value np/n 
and we can calculate a/ap from Eq. (2). If the initial secular eccentricity e < enai, where 
eiim is given by Eq. (16) below, then the eccentricity of the particle is an increasing 
function of time, during the stay of the particle in the exterior mean-motion resonance. 
Eccentricity of the particle can only asymptotically approach to the limiting value eum 
given by the condition X = 0: 

P + <1 _ 1 + 3eiln/2 .,.x 
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Fig. 2. Evolution of eccentricity of dust particle with /? = 0.01 in resonances with a 
planet of mass equal to the Earth's mass, semi-major axis ap — 1 AU and eccentricity 
ep = 0. The left part is for 5/4 exterior mean- motion resonance. Various evolutions 
correspond to different initial values of eccentricity in the resonance. Initial eccentricities 
are: 0.55, 0.45, 0.35, 0.25, 0.15, 0.05. The right part of the figure is for 2/3 interior 
mean-motion resonance. Initial eccentricities are: 0.8, 0.6, 0.4, 0.2. The evolutions are 
numerically calculated from Eq. (4) and are consistent with Eq. (15). 



If the initial eccentricity is greater than eum, then (de/dt) is always negative and the 
eccentricity of the particle is a decreasing function of time, during the stay of the particle 
in the exterior mean-motion resonance. Eccentricity of the particle can only asymptoti- 
cally approach to the limiting value eiim. Characteristic property of the value enm is that 
eiini does not depend on f3. The left part of Fig. 2 depicts evolutions of oscular eccen- 
tricity of dust particle with f3 = 0.01 in exterior resonance 5/4 with a planet of mass 
equal to the Earth mass, semi-major axis ap = I AU and orbital eccentricity ep = 0. 
The evolutions are calculated from numerical solution of Eq. (4) . The first component of 
the right-hand side of Eq. (4) is used as a central acceleration. Asymptotical approach 
to the limiting value of eccentricity eijm ~ 0.2736 (given by Eq. 16) can be easily seen. 
The secular evolution of eccentricity is always a decreasing function of time for interior 
resonances defined by the relation n/np ^ p/{p + q)- Evolutions of oscular eccentricity 
of dust particle with /? = 0.01 in the interior resonance 2/3 {mp = 1 niEarth, ep — 0, 
ap = 1 AU) is shown in the right part of Fig. 2. Evolutions of eccentricities for the given 
resonance and /3 are parallel - the evolutions are shifted along time axis since Eq. (15) 
yields the same value of (de/dt) for the same eccentricity e. li (3 — (e. g., an asteroid) 
we have, from Eq. (15), (de/dt) = 0. 

The crucial question emerges: Are the above presented features typical for real dust 
particles and real physical situations? 

Eq. (15) holds for circular planetary orbit and for the particle with constant optical 
properties in a mean- motion resonance with the planet. Evolution of eccentricity given 
by Eq. (15) will serve as a reference evolution. It will be used as a comparison with 
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the evolutions for the cases of non-circular planetary orbit and non-constant optical 
properties of the particle. 

5. Secular evolution of longitude of pericenter in mean-motion resonances 
with a planet in circular orbit 

Let us assume that a function of the type of Eq. (15) exists for secular evolution of 
the particle's longitude of pericenter if the particle is captured in a mean-motion orbital 
resonance with a planet in circular orbit in the planar case. Let us denote this function 
as S. The assumption is that for a given central star, the planet, the particle and the 
resonance, the function S depends on semi-major axis, eccentricity and longitude of 
pericenter of the particle: 

^)=5(«,e,c.). (17) 

Fig. 3 depicts two evolutions of the orbital elements of the particles with (3 — 0.01 in the 
exterior resonance 5/4 with the planet of mass equal to the Earth's mass, semi-major axis 
ap — 1 AU and eccentricity ep = 0. The particle's initial values of the orbital elements 
for the first evolution, depicted by a solid line, are a w 1.1565 AU (given by Eq. 3), 
e = 0.2, ijj = 90°, and, for the second evolution, depicted by dashed line, a w 1.1565 
AU, e = 0.2, oj = 54°. At the time t = the particles are at the perihelia of their 
orbits. The planet's initial position is at X-axis (axis, from which the longitude of the 
pericenter is measured), in the both integrations. Evolution of the secular eccentricity 
in Fig. 3 is consistent with the behavior expected from Eq. (15). We are interested in 
secular evolution of the longitude of pericenter. Evolutions of the longitude of pcriccnters 
intersect at the time t ~ 614.3 years. Secular values of the semi-major axes are practically 
identical, for the both orbits. The same holds for the eccentricities. However, the values 
of the (doj/dt) differ. This means that the function S{a, e, lu) does not exist for the given 
a, e and lo, since different values of (duj/dt) occur for the same a, e and lo. 

Fig. 4 depicts evolution of the orbital elements of the particle with f3 ~ 0.01 in the 
exterior resonance 5/4 with the planet of mass equal to the Earth's mass, semi-major axis 
ap — I AU and eccentricity ep = 0. Particle's initial values are a « 1.1565 AU, e ~ 0.2, 
uj = 0° (X-axis) for both evolutions. The planetary initial position is at the X-axis in all 
integrations. Various particle's initial positions with respect to the planet yield different 
orbital evolutions. Various values of {dcu/dt) exist during a short time interval at the 
beginning. This means, as in Fig. 3, that function 5* does not exist for the given initial 
values of a, e and u. Although we have chosen only particular values of a, e and lu, we 
conjecture that S{a,e,LLi) does not exist for any a, e and lu. It is impossible to uniquely 
define a function S for various initial conditions r, v, rp, vp yielding a given values of 
a, e and lu. This does not mean that for the given values of a, e and lu there cannot 
exist two various sets of initial conditions leading to the same secular evolution of the 
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Fig. 3. Two orbital evolutions of dust particle with /3 — 0.01 captured in the exterior 
resonance 5/4 with a planet of the Earth's mass, semi-major axis ap ~ \ AU and 
eccentricity ep = 0. Initial values of the first evolution are a « 1.1565 AU, e = 0.2, 
uj = 90° (solid line), for the second evolution a w 1.1565 AU, e — 0.2, ui — 54° (dashed 
line). Evolutions of the longitudes of pericenters intersect at time t w 614.3 years. 

longitude of pericenter. Also in the case (3 — we have found evolutions of the longitude 
of pericenter which have various (cLu/dt) for the same a, e and lu. On the basis of this 
result we conjecture that the function S does not exist even in the case without the P-R 
effect when only perturbation from the planet is taken into account. 



6. Secular evolution of eccentricity and shift of pericenter for noncircular 
planetary orbit 

Fig. 5 shows evolution of semi-major axis, eccentricity and longitude of pericenter of 
spherical dust grain of constant optical properties with /3 = 0.01. The grain is captured 
in the exterior mean-motion resonance 5/4 with the planet of mass mp = 1 rriEarth and 
semi-major axis ap = I AU. Two values of planetary eccentricity are used: ep = and 
ep — 0.2. Black solid line is used for the case ep — 0.2, gray and dashed lines for ep = 0. 
Initial eccentricity of the particle is 0.05. Evolution of the eccentricity for the circular 
planetary orbit is consistent with Eq. (15): evolution asymptotically approaches to the 
limiting value eum ~ 0.2736 given by Eq. (16). Evolution of the particle's eccentricity is 
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Fig. 4. Orbital evolutions of dust particle with (3 — 0.01 captured in the exterior reso- 
nance 5/4 with a planet of the Earth's mass, semi-major axis ap = 1 AU and eccentricity 
ep ~ 0. Initial values of all evolutions are a « 1.1565 AU, e = 0.2, oj = 0°. At the begin- 
ning, the particle is localized in the resonance and its initial position varies with respect 
to the planet. 



a nonmonotonic function of time for the case ep — 0.2. Fig. 5 shows that a limiting value 
of the grain eccentricity may not exist in the exterior resonance, if ep > 0. 

Fig. 6 depicts resonant evolution similar to Fig. 5. The main difference is that the 
black solid line holds for the eccentricity of the planet orbit ep = 0.0167 (eccentricity 
of the Earth orbit), now. Moreover, initial conditions for the particle eccentricity are 
different: the initial value e = 0.2 for the case ep = 0.0167 and e = 0.1845 for the case 
ep = 0. The motivation was to obtain real oscillations of particle's eccentricity around 
the artificial, but analytically solvable (see Eq. 15), case given for ep — 0. Period of the 
oscillations equals to the period of the shift of pericenter (period corresponding to the 
shift in 360°). We have also found oscillatory evolution of the particle's eccentricity for 
an initial value e > eum ~ 0.2736 for the 5/4 resonance. The evolution of the particle's 
eccentricity oscillates around a curve which asymptotically decreases to the limiting value 
eiim ~ 0.2736 (see the left part of Fig. 2). 

If the shift of pericenter is sufficiently fast, then oscillations in the evolution of secular 
eccentricity for exterior resonances exist. Period of this oscillations corresponds to the 
period of the shift of pericenter. This correspondence is caused by libration of conjunc- 
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Fig. 5. Orbital evolution of dust particle with constant optical properties. Particle with 
(3 — 0.01 is captured in the exterior resonance 5/4 with a planet of the Earth's mass, 
semi-major axis ap = 1 AU and eccentricities ep ~ Q and ep == 0.2. Initial eccentricity 
of the particle is 0.05. The case ep — 0.2 corresponds to the black solid line. Gray solid 
line corresponds to the evolution of semi-major axis and dashed line to the evolution of 
eccentricity and longitude of pericenter. The case ep = 0.2 does not lead to asymptotic 
behavior of eccentricity, in contrary to the case ep = 0. 



tions of the planet and the particle around aphelion of the particle's orbit. If the aphelion 
of the particle's orbit is shifted in 360°, then the planetary orbit has the same position 
with respect to the aphelion of the particle orbit (we assume that secular eccentricity of 
the particle does not change significantly). This means that conjunctions of the parti- 
cle and the planet take place in similar trajectories and this leads to similar change of 
the particle's eccentricity. If the shift of pericenter is slow, then oscillations in secular 
eccentricity do not exist. This situation is shown in Fig. 7. 

Fig. 7 compares evolutions of orbital elements of dust particle with j3 ~ 0.01 captured 
in 5/4 exterior resonance with the planet Earth and an "artificial Earth" moving in 
circular orbit. The shift of pericenter in Fig. 7 is much slower than the shift in Fig. 6. As 
a consequence of the slow shift of pericenter, the oscillations in evolution of eccentricity 
are not present in the first 2.5 x 10"'' years for the resonance with the Earth. Evolution of 
the pericenter during the first 2.5 x 10"'' is a nonmonotonous function of time: the initial 



12 



P. Pastor et al.: Motion of dust in mean-motion resonances with planets 




12 3 4 

t[ 10^ years] 




12 3 4 5 

t[ 10' years] 




Fig. 6. Orbital evolution of dust particle with constant optical properties. Particle with 
(3 = 0.01 is captured in the exterior resonance 5/4 with Earth (semi-major axis ap = 1 AU 
and eccentricity ep — 0.0167) - black solid line. Gray color line or dashed line represents 
evolution of particle orbital elements in the resonance with the planet in circular orbit. 
In the case of the resonance with the Earth, the secular evolution of eccentricity exhibits 
oscillations around e{t) given by Eq. (15); period of oscillations corresponds to the period 
of the shift of pericenter. 



decreasing function of time (for approximately 5 x 10* years) is followed by an increasing 
function. 

Our simulations for exterior mean-motion orbital resonances of the first order show 
that maximal values of capture times - for a given (3 and a resonance [j + l)/j - for 
the case ep > are in several tens of percent greater than for the case ep = 0. Also the 
mean capture time for e^ > is greater than the mean capture time for ep = 0. The 
same dependence holds for minimal capture times. This is different from the results for 
nonspherical dust grain (compare Klacka et al. 2005a, 2GG5b; Klacka and Kocifaj 2006a, 
2006b). The greater aspect ratio and smaller volume of the nonspherical grain, the more 
important are nonradial radiation terms (terms with (32 and (3^, in Eq. 41 in Klacka 2004) 
and the capture times in resonances are shorter. If the distance between the central star 
and the planet is greater, then the nonradial radiation terms are less important than 
the effect of the planet - see the term proportional to mp in Eq. (4) - and the capture 
times are greater. These statements for the nonspherical particles are consistent with the 
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Fig. 7. Orbital evolution of dust particle with constant optical properties. Particle with 
/3 = 0.01 is captured in the exterior resonance 5/4 with Earth (semi-major axis ap = 1 
AU and eccentricity ep — 0.0167) - black line. Gray color or dashed curve represents 
evolution of the orbital elements of the particle in the resonance with the planet in circular 
orbit. Initial conditions differ from those used in Fig. 6. The shift of pericenter is much 
slower than the shift in Fig. 6. Evolution of eccentricity exhibits only few oscillations at 
the end of the capture when the shift of pericenter is faster. 



results of Klacka and Kocifaj (2002), Klacka et al. (2005a, 2005b), Klacka and Kocifaj 
(2006a, 2006b). 

As for the interior mean-motion orbital resonances, we have found that noncircular 
planetary orbit can lead to secular increase of the dust grain eccentricity during the 
whole capture time; it seems that ep > ep{critical) > is required, where ep(critical) 
depends on the type of the resonance, /3 and mass of the planet. This is different from 
the case ep == 0, when only secular decrease of eccentricity exists (see Eq. 15). The 
situation is illustrated in Fig. 8, which holds for resonance 2/3, f3 — 0.01, ap = 1 AU, 
ep = 0.2 and planetary mass equal to the Earth mass. The obtained result is similar to 
the evolution of nonspherical particle with circular planetary orbit (compare Fig. 8 with 
Fig. 2 in Klacka et al. 2005b, or, with Fig. 3 in Klacka and Kocifaj 2006b). Fig. 8 shows 
also very complicated behavior of the shift of pericenter of the particle. Our simulations 
show that capture times of particles in the interior resonances are larger for ep 7^ 
than for the case ep — 0, if initial particles eccentricities are small. When eccentricity 
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Fig. 8. Secular evolution of orbital elements of dust particle for the interior resonance 
2/3 with a planet. Spherical particle of constant optical properties is characterized by 
(3 — 0.01. Mass of the planet equals the Earth mass, semi-major axis ap — \ AU and 
planetary eccentricity ep — 0.2. Secular evolution of particle eccentricity is an increasing 
function of time during the capture in the resonance. Shift of pericenter is similar to the 
case ep — 0. 

of the particle in an interior resonance with the planet in circular orbit decreases to 0, 
the capture is ending because of lack of positive orbital energy given by the planet (Liou 
and Zook 1997). However, the eccentricity of the particle can be an increasing function 
of time, if the particle is captured in an interior resonance with the planet moving in 
elliptical orbit. The particle can approach to the planet and gain a sufficient amount of 
positive orbital energy to prevent a decrease of semi-major axis caused by the P-R effect. 
This is explanation of the longer capture times for ep 7^ 0, if initial particle eccentricities 
are small. 



7. Influence of changing optical properties on evolution of orbital elements 
in mean-motion resonances 

We are interested in orbital evolution of a particle with changing optical properties. As 
a central acceleration we will use —GMq{1 — f3o)ep/r'^, where /3o is the value of f3 at the 
time t — 0. 
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We will calculate time derivative of semi-major axis and eccentricity from perturba- 
tion equations of celestial mechanics for Eq. (4) without perturbation of a planet. We 
get 

dt 1-el \^GMQ{l-f3o) H 



-/3o 

2ap„ I3GMq 



1 — el cr'^ 

pa 



[2(e/3oSin/^„)2^(l + e^„cos//3„)2] ^ (ig) 



- (/?-/3o) — J— sm//3o 



dt \IGMq{1-Po) 

^^ [e/3o sin^ fpo + 2(eft, + cos fp^)] , (19) 

where ap^, ep^ and fp^ arc oscular semi-major axis, eccentricity and true anomaly, pp^ = 
Q/Sg (l — ej) ) is parameter for elliptical orbit and r = p^|-,/(l+e^„ cos fi3„). Secular evolution 
of semi-major axis and eccentricity can be obtained by averaging over one orbital period 
of the type 

(5) ^^ / \'g{fpM0o , (20) 

where g is any quantity. (3 is even function of true anomaly fp^. Thus, the average value 
of the first terms in Eqs. (18)-(19) is zero. We can write 

/dai3o\ ^ 1 GMq 

\ dt / ^a^„(l-e|J3/2 , 

X / /? [2{e0„ sin fp,f + (1 + e^, cos fp,f] dfp, , (21) 



def3„ \ _ 1 GMq 

— X 



X / P [ep„sm^ fp„+ 2(ep„+ COS fp„)]dfp„ . (22) 

JO 

P is approximately constant in Eqs. (21)-(22), see Fig. 1 (for more details about this 
approximation see Klacka et al. 2007). We get for secular evolution 

da«„ \ „ GMa, 2 + 3e2 



d(^0o\ _ 5 GMq ep„ 

Po — — -o-n 9 M/9 • i-^V 



dt 2"^" c al (l-e2 )i/2 

Eqs. (23)-(24) are identical to Eqs. (13)-(14). Inserting Eqs. (23)-(24) into Eq. (12) (with 
the assumption that Eqs. (23)-(24) hold also for period of resonant oscillation of semi- 
major axis), we obtain equation which is identical to Eq. (15). This means that the change 
of optical properties does not significantly influence the evolution of secular eccentricity 
in mean-motion orbital resonances. 
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Now, we will derive an expression for secular evolution of the longitude of pericenter 
for Eq. (4) without action of a planet. Perturbation equations of celestial mechanics yield 

- (P - Po) ^r~ COS //3o + 



dt efi, "^ GMq{1-Po) 

1 GMq ^ ^ , ^ 

+ /3 ^ sm /^„ e/3„ cos //3o - 2) . 25 

After averaging, using also Eq. (20), we have 
dujpo \ _ GMq I pjr„ 



X / P COS fpjjp^ . (26) 

Jo 

Since (3 is an even function of true anomaly /^„ , the second term in the right-hand side 

of Eq. (25) equals zero, after averaging. Also the average value of the term proportional 

to /3o is zero. The shift of pericenter caused by the change of optical properties is in 

the same direction as the particle orbits the Sun, since (3 is an increasing function of 

heliocentric distance (see Fig. 1). 

Fig. 9 depicts secular evolution of orbital elements of the particle with changing 
optical properties. Solid black line is used for spherical particle with radius i? = 5 /im 
and density p = 2 g/cm'^ captured in 4/3 exterior resonance with the Earth. Gray line 
or dashed line is used for evolutions of the particle captured in the resonance with the 
planet of mass equal to the Earth mass, semi-major axis ap — 1 AU and eccentricity 
ep = 0. Evolution of eccentricity in Fig. 9 is similar to evolution of eccentricity in Fig. 6. 
Evolution of the eccentricity for the case ep — Q asymptotically approaches to the limiting 
value eiini ~ 0.3108 given by Eq. (16). Evolution of eccentricity for the resonance with the 
Earth exhibits oscillations around the evolution for the resonance with circular planetary 
orbit. Period of oscillations corresponds to the period of the shift of pericenter. Evolution 
of the longitude of pericenter exhibits an increase in time. This is the difference from the 
behavior presented in Fig. 6. Behavior of the longitude of pericenter in Fig. 9 is caused, 
probably, by the influence of the changing optical properties, since the motion of the 
particle is prograde. 

Fig. 10 compares evolution of a particle of constant optical properties with the evo- 
lution of another particle of changing optical properties. Both particles are characterized 
by the same initial conditions, as for the orbital elements and position with respect to the 
planet. Mass of the planet is equal to the mass of the Earth, semi-major axis ap = 1 AU 
and eccentricity ep = 0. The particles are captured in 3/2 exterior resonance with the 
planet. Fig. 10 presents that the shift of pericenter is in positive direction/orientation, 
for the case of changing optical properties of the particle. This is in coincidence with 
Eq. (26), since the motion is prograde. The shift of pericenter is in negative direction for 
particle with constant optical properties. 
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Fig. 9. Orbital evolution of dust particle with optical properties varying with heliocentric 
distance; radius of the particle R = b /ini, density p ~ 2 g/cm^. The particle is captured 
in the exterior resonance 4/3 with Earth (semi-major axis ap = 1 AU and eccentricity 
ep — 0.0167) - solid black line. Gray color or dashed curve represents orbital evolution 
of the particle in the resonance with the planet in circular orbit (ep = 0). In the case 
of the resonance with the Earth, the evolution of eccentricity shows oscillations around 
e{t) given by Eq. (15), with period of the shift of pericenter. Greater part of the shift 
of pericenter exhibits an increase in time and this differs from the behavior presented in 
Fig. 6. 



It is also possible to find the shift of pericenter in negative direction for particle with 
radius R = 2 /im and density p = 2 g/cni^. This case is depicted in Fig. 11. However, 
the shift of pericenter in negative direction is rare, for this type particle: the situation 
happens for initial conditions shown in Fig. 12. 

As for the shift of pericenter, we may use perturbation equations of celestial mechan- 
ics. Central acceleration will be —GMq{1 — Po)sr/'>''^, where f3o is the value of /3 at time 
t = 0. Subtracting the central acceleration from the right-hand side of Eq. (4) yields the 
perturbation acceleration. Again, the planar circular restricted three-body problem and 
the P-R effect are considered. We get 



1 



P/3o 



dt ep, VGMo(l-/3o) 

GMq 1 1 



GMq 

(/3-/3o) — 5— cos ff3. 



P _pi 1. ^ (e^ cos //3„ - 2) sin /^„ 
r2 c e«„ 
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Fig. 10. Black line depict orbital evolution of dust particle with optical properties de- 
pendent on heliocentric distance; radius of the particle R ~ 2 /im, density p = 2 g/cm'^. 
Particle is captured in the exterior resonance 3/2 with a planet of Earth mass, semi- major 
axis ap = 1 AU and eccentricity ep — 0. The shift of pericenter is in positive direction. 
Gray color or dashed curve represents orbital evolution of the particle in the resonance 
with the planet (identical initial conditions), but optical properties of the particle are 
constant. The shift of pericenter is in negative direction for particle with constant optical 
properties. 



a/3o(l 



Y = - 
Z 



Gmp 

~^ V GMq (1 - /3o) 

cos //3o 



^0O> ^y 



\r — rp\ 
sin //3o 



rp 



|r-rp|3 r?, 



X Z 



■ sin (6^30 - 0p) + cos (Op - ujfj^) 



(27) 



1 + e/3o cos //3o 

where Of3g — ajp^ + //3„ is the position angle of the particle on the orbit (measured from 
an X-axis) and Op is the position angle of the planet on the orbit (measured from the 
X-axis) . 

Initial conditions depicted in Fig. 12 are: /3 = /3o, //Jq = 0°, 0/3^ = 0°, uji3„ = 0°, 
9p = 0°, and r is a little less than rp. Inserting the initial conditions into Eq. (24), one 
obtains duipg/dt < °/year. If r is a little greater than rp, then dui/^g/dt > °/year. 
When the particle, is approximately, in the position depicted in Fig. 12, then the first 
term in Eq. (27) is negative - /3 is an increasing function of heliocentric distance and 
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Fig. 11. Evolution of orbital elements of dust particle which optical properties depend 
on heliocentric distance. Radius of the particle R ~ 2 /im, density p = 2 g/cm"^. The 
particle is captured in the exterior resonance 2/1 with a planet of mass equal to the 
Earth mass, semi-major axis ap — 1 AU and eccentricity ep ~ 0. Shift of pericenter is in 
negative direction for t < 400 years. This situation happens for initial conditions shown 
in Fig. 12. 



Po is the value of (3 at perihelion of the particle's orbit (i = 0). The sign of the second 
term depends on the sign of sin/^g and the third term is negative when the particle's 
position corresponds to that depicted in Fig. 12. The third term is dominant and the 
shift of pericenter is negative, since /3 « /3o, ff3„ ~ 0° and the particle is near the planet. 
The oscular evolution of the longitude of pericenter is negative when the particle is in the 
interior part of the planetary orbit and the conjunction is near the particle's perihelion. 
This effect is relevant also for the secular evolution of the longitude of pericenter (shift of 
pericenter): when the particle moves outside the planetary orbit, then oscular evolution of 
the longitude of pericenter may be positive, but averaging over the orbital period yields 
negative sign. This may be explanation of Fig. 11 which depicts a decreasing secular 
evolution of the longitude of pericenter. 

If the motion of dust particle is retrograde and the motion of the planet is prograde, 
then the shift of the particle pericenter should be, mainly, in negative direction, for the 
changing optical properties of the particle (Klacka et al. 2007). However, this depends 
on the detail behavior of the change of optical properties. 
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Fig. 12. Initial conditions for orbital evolution shown in Fig. 11 in the orbital plane 
XY. Ellipse represents trajectory of dust particle during the first two years. Sun is at 
the origin. Circle represents trajectory of the planet. Black cross represents the initial 
conditions of the particle. Black dot on the right represents initial position of the planet. 
The particle is initially inside the planetary orbit, but it is captured in the exterior 
resonance. 



8. Conclusions 

The contribution deals with the effect of solar electromagnetic radiation on dynamics 
of cosmic dust particles in mean-motion orbital resonances with a planet of mass equal 
to the mass of the Earth. We discuss not only the planar circular restricted three-body 
problem and the Poynting- Robertson effect with constant optical properties of the spher- 
ical particle. We admit also nonzero eccentricity of the planet and the change of particle 
optical properties with heliocentric distance. The paper concentrates on important prop- 
erties of motion of dust grain in the zone of mean-motion resonances. We are interested 
mainly in pericenter motion and evolution of eccentricity in the resonances. 

Our numerical integrations suggest that any analytic expression for secular time 
derivative of the particle's longitude of pericenter does not exist, if a dependence only 
on semi-major axis, eccentricity and longitude of pericenter is considered (P-R effect and 
mean-motion resonance with planet in circular orbit are taken into account). 

If planetary eccentricity is close to zero and the shift of pericenter is sufficiently fast, 
then oscillations of dust grain secular eccentricity exist. The oscillations occur around the 
curve corresponding to secular evolution of the grain eccentricity in the planar circular 
restricted three-body problem with the P-R effect. This holds for exterior mean- motion 
orbital resonances. Nothing like this was found in the case of interior resonances. However, 
interior resonances can exhibit systematic increase of secular eccentricity of the grain 



p. Pastor et al.: Motion of dust in mean-motion resonances with planets 21 

during the capture in the resonances. This is true when eccentricity of the planet is larger 
than some critical value depending on the type of the resonance, mass of the planet and 
optical parameter (3. The case ep = yields only secular decrease of eccentricity, for 
interior resonances. 

Our numerical simulations show that noncircularity of the planetary orbit stabilizes 
motion of spherical dust grain in the mean-motion orbital resonances. Maximal capture 
time (and also mean capture time for many captures in numerical simulations) of the 
grain for a given resonance and particle is greater than it is for the case of circular 
planetary orbit. This holds both for exterior and interior resonances. Nonsphericity of 
the grain destabilizes motion in the resonances: the greater aspect ratio and smaller 
volume of the grain, the shorter capture time. 

If a change of optical properties of the spherical grain with heliocentric distance is 
also considered, then the shift of pericenter is dominated in positive direction/orientation 
for prograde motion of the particle; this holds both for circular and noncircular planetary 
orbits and exterior mean- motion orbital resonances (see, e. g.. Fig. 10). If planetary orbit 
is characterized by large eccentricity, then secular evolution of dust grain eccentricity 
may exhibit complicated behavior. 

Spherical dust grain in the planar circular restricted three-body problem with the 
Poynting- Robertson effect is characterized by a monotonic secular evolution of the grain 
eccentricity (this is true both for constant and variable optical properties of the grain) , 
in the exterior mean-motion orbital resonances. Such kind of behavior does not exist if at 
least one of the above mentioned assumptions - sphericity of the grain or circular orbit of 
the planet - is cancelled. As for the interior resonances, the circular planetary orbit yields 
secular decrease of the eccentricity, while nonsphericity of the grain or noncircularity of 
the planetary orbit may yield secular increase of the eccentricity, also. There is some kind 
of unification of qualitative kinematical behavior of interplanetary dust grains under the 
action of more real physical forces: there is not great difference between evolution of 
orbital elements of spherical and nonspherical dust grains if more general physical forces 
are taken into account. Our results show that gravitational effect can mimic nongravita- 
tional effect. However, the following question is still unanswered: Can spherical grain be 
captured into a resonance when the secular evolution of particle semi-major axis is an 
increasing function of time? 
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